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METHOD AND APPARATUS FOR PREDICTING PLASMA-PROCESS 

SURFACE PROFILES 

FIELD OF THE INVENTION 

TUsimrendon relates to plasm processing of smiconductordev^ In 
particular, this invendon provides method and apparatus for predicting the sui&ce profile 
that a ^en plasma process will create. 

BACKGROUND OF THE INVENTION 

The manufacture of a semiconductor device routinely entails a sequence of 
processes in each of which a substrate is CTcposed to a partially ionized gaseous medium in 
order to effect deposition of material onto the substrate or removal of material therefrom. 
Such a plasma process is defined by an ensend}le of macroscopic hiput 
parameters — including, for example, pow^, temperature, pressure, inlet gas composition, 
substrate material, duration of exposure— eveiy one of which affects the resulting profile 
on the substrate. Traditionally, a suite of values of these input parameters suitable for 
creating a given set of device features has been determined by trial and error. 
Devdopment of a single process by tins enqnrical approach is costly and time-consuming, 
requiring treatment of several patterned wafers and subsequent study of the resultix^ 
profiles by scanmt^g electron microscopy. Because of the unpredictable way a small 
change in one input parameter may affect the profile, any modification of the layout — ior 
example in device dimension, pattern dend^ on the wafer, diange in total open 
area— fix>m one implication to another has often necessitated redevelopment of the 
process, with the attendant outiay of resources. 

Recent advances in device &brication tedmology are rendering this q)proach even 
more onerous. Decreasing feature sizes demand tighter tolerances on feature dimensions 
and morphologies, so that the number of trials required to optinuze a given process is 
inoeasing The acceleration of wafer diameter growth and the con^lete redesign of the 
process involved with an incremental change in diam^er have increased the xumiber of 
times tills empirical process must be repeated. The increasing use of devices tailor-made 
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to a specific application also increases the amount of development and optimization 
activity required. 

An altOTiative, computational approach would derive input parameters from a 
complete physical description of a plasma process inducing a plasma model for describing 
the coupling between the macroscopic input parameters and the maCTOScopic fluxes, 
concentrations and enei^ distributions of the various species in the plasnu^ and a profile 
amulator for atomistically deteramiing from the nuu^oscopic fluxes the resulting etch or 
depoation rate along the wafer surface and calculating the profile evohition therefix>m. 
Ideally, such a physical description of plasma etching and deposition processes would 
enable the ab initio selection of the maax)scopic mp\jx parameters appropriate fi>r 
generating a desired profile on the substrate, eliminating the need for expenave and time- 
consuming test sequences. 

Research in tlus field has done much to duddate mechanisms at work in plasma 
processes, and thus has contributed scaling laws that could frame a physical descriptioa 
However, notwithstandmg the availability of computational means suffidentiy powerfiil to 
perform the necessary calculations based on known scaling laws, the in^ylemcntation of 
such an ab initio approach has been limited by lack of data. For example, the mann^ in 
which the vahies of some coeffidents in these laws depend on the particulars of a gwen 
process has been imknown as yet In past investigations, determination of the value of 
such a scaling coeflBdent consistent with a plasma process defined by a givra set of input 
paramet^s has typically been done by comparing a fiiushed profile, created by applyirig 
that process, with a simulated profile mcluding one or more of these coeffidents as 
adjustable parameters. Such Wndright evaluation may promote understanding a ^en 
coeffident*s role m a scaling law, but it has not afforded the ability to predict profile 
evohrtion for any process defined by a set of izq>ut parameters diflFering from the set used 
in tiie experimental process used to derive the value of that coefficient. 

SUMMARY OF THE INVENTION 

The invention enables prediction of a process smface profile that a plasma process 
of interest, defined by a set of process values of input parameters including input 
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variables, wiU create on a process substrate. The prediction is based on a test surface 
profile, defined to be the experimental outcome of a test process defined by a set of test 
values of input parameters, mcluding the input variables, such that at least one of the input 
parameter vahies defining the test process differs from the corresponding input parameter 
value defining the process of interest. (As used herein, the term "input variable" 
0 designates an input parameter whose vahie is free to vary within the context of a giv^ 

application of the invention to a plasma process.) In another aspect, the invention 
provides a technique for determining, using a test surfece profile, process values of input 
variables that wiU produce a desfred surfiice profile. In this aspect, the test surfiice profile 
differs fiom the desired profile, or the value of some fixed input parameter defining the 
,5 test process differs from the corresponding parameter value defining the process for 
producing the deared profile. 

In related aspects, the invention also provides apparatus for predicting a process 
8uifiM» profile and determiniiig process vahies, a method of configuring a plasma reactor, 
a method of making semiconductor devices requiring limited empirical calibration, and a 
10 device made according to the method. 

The semi-empirical approach of the invention incorporates a caBbration procedure 
that operates as follows. An initial mathematical surface profile model, formulated in 
terms of fixed input parameters, input variables and a phirality of coefficients of unknown 
values, is evaluated at the test vahies of the input variables to generate a quantitative 
j3 approximate prediction of the profile created by the test process. The vahies of the 
unknown coefficients are then optimized to minimize the difference, as determined by a 
quantitative comparison, between tiie test suifiwe profile and approximate profile 
prediction. By substitution of tiiese optimum values for tiie unknown coefficients in the 
initial surfiice profile model, the cafibration provides a final matiiematical model having 

30 the input variables as unknowns. 

The calibrated final matiiematical modd can be used in any one of several aspects 

of the invention. Inserting process values of interest into the final mathematical model 
enables prediction of tiie process surfece profile tiiat would be created by applying to a 
process substrate tiie plasma process sequence defined by tiie process vahies. The 
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process values used for this predictive function may be constant, corresponding to a single 
static process or may vary over tirae> for example representing a multistep process. 
Ahematively, process values appropriate for creating a desired suifece profile are derived 
from the final mathematical profile model and the desired profile. 

the calibrated final mathematical model is also usefiil in formulating a new plasma 
process sequence using a different reactor de^ than that used to a*eate an available test 
sur&ce profile, or having one or more process vahies sufiBciently different from the 
corresponding test value to place the new sequence in a regime for which the test process 
has limited predictive power, in such cases, input variable process values derived using the 
calibrated coefficients may s^e as a sensible starting point for empirical optimizafion, so 
that the number of wafers that must be spent in finding the ultimate process sequence is 
reduced. 

Thus by using data available from direct measurement of a limited number of 
laboratory-generated mscrostructures, the invention provides powerfiil techniques for 
plasma process device manu&cture that incorporate agnificantly less experimentation 
than fiiUy empirical methods, yet do not require the body of fiindamental data needed for 
a complete computational simulation. 

The semi-empirical approadi of the invention is not limited to any particular 
macroscopic physical description of the sur&ce profile or underi^g pla^a modd or 
profile simulator. The mathematical models may incorporate scaling laws derived from 
basic plasma physics and chenustiy, numerical emulations or expoiment. In general, the 
plasma model is configured to return fluxes^ energy and angular distributions arriving at 
the substrate surface fi>r spedes of mterest, based on parameters related to the fimction of 
the reactor. The plasma model may entail a conq)lete plasma simulation tool with all 
known cross sections for the gas-phase chemistry and plasma sur&ce interactions, or 
merely a look-up table of experimentally measured fluxes and observed trends. In either 
approach, the model graerally includes at least one coeffident of unknown value that is 
determined by the calibration of the inventioa Individual coeffidents of unknown value 
in the plasma model provide reference magnitudes or characterize the coupling between 
input parameters and aspects of plasma behavior, and thus they collectively describe the 
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5 physical properties of the medium generated by the plasma process. Alternately, the 
plasma model may be eliminated altogether as a distina module, the profile simulator 
treating the fluxes as adjustable coeflBcients of unknown value to be determined by the 
calibration. 

The profile emulator preferably, for accuracy and versatility, calculates the 
10 movementof each point of the substrate sur&ce. An appropriate profile sumilator in 
general includes a local transport module for calculating local fluxes arriving at ead\ point 
of the substrate sur&ce, a ate balance model for calculating local rates of etching and 
depo^on, and a surface advancement algorithm for translating these mechanisms into net 
surface movement. In a preferred embodiment, the profile simulator derives the local 
15 fluxes from the macroscopic parameters provided by the plasma model using Monte Carlo 
methods; a Langmuir-type model describes the ate balance; and the surface advancement 
algorithm incorporates a shodc-firont-traddng approach. Generally speaking, each of 
these three segments contributes coefiBdents of unknown value to be determined at 
calibration. The coeflBcients of iinknown value in the profile shnulator characterize 
20 couplings between the processing medium and the behavior of the substrate surfece and 
thus collectively describe the substrate sur&ce evolution. 

In a preferred embodiment, the quantitative comparison of the calibration 
compares the test surface profile wth the approwmate prediction at each point along the 
substrate sur&ce, and the local difiPerences constitute a residual fimctioa A 
23 multidimenaonal nonlinear least-squares technique, sudi as a modified I^enberg- 

Maiquandt algorithm, is used to determine how each of the unknown coefiBdents aflFects 
the residual fiinction and then to adjust the unknown coefficioits so as to attempt to 
maximize some figure of merit and provide acceptable goodness of fit. 

As used herein, tiie tenn **test profile" denotes all of the experimental profile data 
30 available for a given calibration and may include multiple snapshots, each taken after a 
diflferent length of procesang time but otherwise under the same process conditions as the 
other snapshots in the profile. Or, more generally, the snapshots in a given test profile 
may be generated using test processes var^g in aspects other than procesang time. 
Also, this desaription of the invention emphasizes the general case of a plasma process 
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defined by a phirality of input variables and mathematical surface profile models including 
a plurality of coefficients of unknown value. However, the invention also encompasses 
applications to treatment sequences having only one free input variable and/or 
mathematical models for \dvich the value of only one coefficient is unknown. 

BRIEF DESOOPTION OF THE DRAWINGS 

The invention description below refers to the accontqpanying drawiqgs, of i^ch: 

FIG. 1 is a flow chart illustrating the calibration sequence for the method of the 
invention; 

FIG. 2 schematically depicts sur&ce kinetics according to the Langmuir model ; 

and 

FIG. 3 schematically illustrates a representative hardware environment for the 
present invention. 

DETAILED DESCRIPTION OF AN ILLUSTRATIVE EMBODIMENT 

The general procedure for performing the calibration of the invention based on a 
^ven empirical test process is illustrated in FIG. 1. In a first step 200, test values of the 
input parameters are selected. The mput parameters include descriptors of the pre- 
processing state of the substrate ("substrate parameters") and quantities definuig the 
operation of a giv^ plasma reactor during the test process ("reaction parameters"). The 
substrate paramet^s may include, for example, gross dimmrions of the substrate, the 
distribution and dimensions of any features on the sur&ce, and the substrate compo^tion. 
The reaction parameters may inchide any of power levels, temperature, pressure, inlet gas 
compoation and length of time the substrate is treated by the plasma process. 

The test values of the reaction variables and fixed reaction parameters are used m 
step 210, in which a plasma model, containing a mathematical description of the plasma in 
terms of the reaction parameters for the reactor of interest and coefficients of unknown 
value, operates to approximately determine the macroscopic properties of the plasma 
generated by the test process. The approximation carried out in step 210 is based on 
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rough preliminary values of the plasma model coeflBcients of unknown value, provided in 
stq)215. 

In step 220 a profile simulator, containing a mathematical model of the time 
evolution of the substrate surfece profile in terms of the substrate parameters, the 
macroscopic plasma properties and coefEdents of unknown value, uses the results of step 
210 and the substrate parameters provided in step 200 to approximately predict the effect 
of the test process on the substrate sur&ce profile. The approximation carried out in step 
220 is based on rough preliminary values of the profile amulator coefficients of unknown 
value, provided in step 22S, Thus, the cumulative output of steps 200 through 220, 
particularly of the operation of the plasma model in step 210 and the profile simulator m 
step 220, which together comprise an mitial mathematical sur&ce profile model, is a 
quantitative but approximate prediction of the profile created by the test process. 

In step 230 the test values of the input parameters provided in step 200 are used 
to provide a test surface profile, created experimratally by subjecting a test substrate to a 
test process, in the reactor of interest, defined by the test values of the input parametm. 
This test sur&ce profile is quantitatively compared in stq> 240 with the quantitative 
approximate prediction resulting firom step 220. The difference between the test sm&oe 
profile and the appro^dmate prediction is evaluated according to some criterion applied in 
step 240. In general, on the first pass, in whidi the approximate prediction is computed 
using preliminary values of unknown coefiBdents, the residual is not suffidentiy small to 
satisfy the criterion, and the calibration procedure advances to step 250, in vMd^ the 
values of the unknown coeffidents are adjusted so as to mitigate the re^dual. The 
adjusted values are then passed to the plasma model and/or the profile simulator throu^ 
st^s 215 and/or 225, for anotiier iteration of tiie steps that render the approximate profile 
prediction and of the comparison in step 240. Iteration continues until the discrepancy 
between the test surface profile and the approximate prediction is adequately minimized. 

Thus the iterative comparison of step 240 and coeffident adjustment of step 250 
effect a calculation of optiiinim values of the coefiBdents of unknown value appearing in 
the initial mathematical surface profile modd. After the calibration of the invention 
deteniunes the coeffident values suitable for the mathematical sur&ce profile model, the 
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resulting final mathematical model can be used to account for the effect on a substrate of 
plasma process sequences differing fi-om the test process in one or more process variable 
values, so long as the semi-empirical model properly accoimts for factors such as loading 
and aspect ratio dependencies.- 

Preferably, the approximate profile prediction generated by the profile simulator 
comprises a series of frames, each corresponding to a regular tune increment^ in step 220 
and the compare step 230 only compares frames of the approximate profile prediction that 
correspond to cumulative exposure times roughly equal to the time of the snapshots in the 
test sur&ce profile. If the test profile includes multiple snapshots at different esq^osure 
times and/or at a different test vahie of one or more input parameters other than time, step 
240 compares each test snapshot with the appropriate firame of the profile prediction and 
the system operates to minimize the residual over the entire pairwise comparison. 

As will be evident to those skilled in the art, many configurations dq>arting firom 
the procedure shown in FIG. 1 fall within the scope of the inventioit For example, the 
entire initial mathematical surface profile model may re^de m a aoi^e module rather than 
bmg divided into distinct parts 210 and 220 to deal separatdy with phenomena acting 
over different length scales. Or, depending on the nature of the initial mathematical 
surface profile model and con:q>arison algorithm used, insertion of the test values and 
preliminary coeffidents into the mitial mathematical surface profile model may be delayed 
until the comparison step. 

In a preferred embodiment, the plasma model characterizes the spedes of interest 
only as belonging to a general class of actors, for example, as bdng either charged 
particles, and thus attracted to the substrate by any applied bias, or neutral species sudi as 
gas molecule and exdted radicals. Physical modds--for example. Maxwell's and 
Boltzmann's equations— underiying the functional dependence m Xems of goveniing 
experimental input parameters of such plasma descriptors as partide fluxes, energy, and 
angular cUstributions are well imown in the art (See, e.g.» lieberman and Licfatenberg, 
Principles of Plasma Discharges and Materials Processing, John Wiley, 1994.) Much is 
known m the art about how these descriptors scale with the input parameters, based on 
such basic physical nk)dds in conjunction vrith experimental data. However, the absolute 
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5 value of these fluxes or distributions is not known a priori for a given process. Any such 
coefficient in the model whose value is unknown is treated by the calibration as an 
adjustable parameter. 

In general, the profile simulator includes a local transport model for calculating 
fluxes arriwig at each point along the substrate surface, a fflte balance model for 
10 calculating the resulting local etdi and dq>osition rates, and a surfece advancement 
algorithm. In a preferred embodiment, a Monte Carlo approach is used for tiie local 
transport, for which information from atomic scale, rinmlation is appropriate. 

The site balance segment is preferably based on a Langrauir-type model of the 
kinetics of particle-surfece interactions. The model characterizes spedes that remove 
13 material from the substrate as etchants and those tiiat deposit material onto the substrate 
surface as inhiT)itors. With reference to HG. 2, r, and respective^ 
of etchants and inhibitors arriving at the surfitceofthe substrate 270. Both etchants and 
inhibitors react with the surfece when they encounter a vacant ate. Etdiant and inhibitor 
spedes react at the surface wfth respective open-fidd reaction coeffidents of and 
20 to form etchant site 275 and inHbitor dte 280. The overall occupandes of the substrate 
surface sites resulting from reactions witii etchant and inhibitor spedes are represented 
respectively by fractions 0, and ^^,whhtiie number of vacant dtes 285 bdiig equal to 
1^0^^$^^ Etdi products can form and thOTiallydesort) from the surface, lonsarriwig 
at the surface, at flux , can remove etchants tiiat have reacted with the surfece by 
23 physical sputtering and ion-enhanced etching; arriving ions can also remove inhibitors 
from the surfece by physical sputteriiig. Taking into account these assumptions and 
mechanisms, the site balance, in the steady state, tor tiie two types of occupied sites be 
expressed respectivdy as 

a.^^T,S,(l-0,-9yeX,Tu-W^-ka,A^O and 

at 

wherein , , and k^^ are coeflidents, generally of unknown value, assodated witii 
the ion-enhanced, phyacal sputtering and thermal etching mechanisms, respectively. The 
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parameters and are products of ion-enhanced etching and physical sputtering 
yields, respectively, with ion flux, integrated over incident ion ener^es greater than 
respective threshold energies, expressed respectively as Eu^j and JE^^ , and over all 

angles, 

llY^{^,E)T,{^,E)d^= j dE{E^-EoS\id4W>EW'E)- 

We have assumed that the etch ^dds are products of angular fiincdons and the 
square root of ion energy, which dependence has been observed experimentaUy. 
However, other scaling laws could be used instead. The yield fiinctions only represent 
fimctional dependences, the absolute magnitudes bdng lumped into the A:'s along with 
other constants. 

Expressions for 0, and 9^ can be derived from the steady-state site balance 
equations on each type of surface present on the substrate, for example the substrate top 
and trench bottoms and sidewaUs. At each point on the sur&ce, the etdi rate can be 
written 

ER=k^ r,+^A2r,,+*^F.. 

The coeffici^ts y ^n* K y^^^ constants assodated with physical 

sputtering, ion-enhanced etdung, thermal etching and ion-induced deposition, 
respectively. Incorporating the ^ressions for 0^,0^^ F^ and F^ renders the etch rate 
in terms of the plasma characteristics and the coefiBdents k , each of whidi is in prindple 
an adjustable paratiieter that can be deternuiied by the calibratio 

It is also preferred to employ an analytic scheme for surface advancement so that 
the fine features can be resolved more accurately. One such sdiemes known in the art and 
capable of modeling fine feature aspects, such as sharp comers, is the method of 
characteristics, also known as shock-front-tracking algorithm. (See, e.g., S. Hamaguchi, 
''Modeling of Fibn Deposition for Kficrodectronic Applications'", JJan Films, vol. 22, p. 
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5 81, S. Rossnagel, ed.. Academic Press, San Diego, 1996). Another is the level set 

approach. (See, e.g., J. A, SetWan, Level Set Methods : Evolving Interfaces in Geometry, 
Fluid Mechanics, Computer Vision, and Materials Science, Cambridge University Press 
1996.) The shodc-front-tracking approach models the surface OLe., the interfece layer 
between vacuum and solid) as a collection of piecewise continuous line segments^ for 

10 each ofwhich a rate ofmovement is calculated. Thepossibility of each segment's 
advancing or receding along its normal independently of the movement of the other 
segments aUows for multiple potential solutions for the resulting sur&ce. In order to 
avoid multiple sohitions, these analytic schemes modd the points between the line 
segments as shocks (te,, discontinuities in the slope), and properly track the movement of 

15 the shodcs. 

Refer now to FIG. 3, which iDustrates, in block-diagram form, a hardware system 
incorporating the invention. As indicated therein, the system includes a system bus 355, 
over which all system components communicate, a mass storage device (such as a hard 
disk or optical storage unit) 357 as well as a main system memory 360. 

20 The operation of the illustrated system is directed by a central-processing unit 

("CPU") 370. The user interacts vwtii the system using a keyboard 380 and a position- 
sensing device (e.g., a mouse) 382. The output of eitiier device can be used to deagnate 
information or select particular areas of a screen display 384 to direct fiinctions to be 
performed by the systm. 

23 The main memory 360 contams a group of modules that control the operation of 

CPU 370 and its interaction with the other hardware conq?onents. An operating system 
390 directs the execution of low-level, basic system fimctions such as memory allocation, 
file managonent and operation of mass storage devices 357. At a Mgher level, an analysis 
module 392, implemented as a series of stored instructions, directs execution of the 

30 primary functions performed by the invention, as discussed below. Instructions defining a 
user interface 394 allow straightforward interaction over screen display 384. User 
interfiice 394 generates words or gn^)hical images on display 384 to prompt action by the 
user, and accepts user commands firom keyboard 380 and/or po^on-sensing device 382. 
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The main memory 360 also includes one or more databases 396, in general cont^iung any 
of the test or process values of input parameters including input variables, the desired 
profile, the test surface profile and rough preliminary values of the coefficients of 
unknown value in the plasma model and profile simulator. 

It must be understood that although the modules of main memory 360 have been 
described separately, this is for clarity of presentation only; so long as the system 
performs all necessary fimctions, it is inunaterial how they are distributed within the 
system and its programming architecture. 

The test sur&ce profile is produced experimentally, as is weO known in the art, by 
subjecting one or more test substrates to a test process in a plasma reactor and measuring 
the resulting surface profile using, for example, scanning electron microscopy. The 
desired and test surface profiles may be supplied to the hardware system in electronic 
format or as graphic hardcopy, in which case the image(s) is processed by a digitizer 398 
before rmmerical comparison with the approximate prediction. The digitized profile is 
sent as bitstreams on the bus 355 to a database 396 of the main memory 360. The test 
surface profile may be stored in the mass storage device 357 as well as in databases 396. 

As noted above, execution of the key tasks associated with the present invaition is 
directed by analy^s module 392, vMch governs the operation of CPU 370 and controls 
its interaction with main memory 360 in performmg the steps necessary to provide a final 
mathematical surface profile model including calibrated optimum values of the unknown 
coeffidents m the initial sur£tce profile model; and, by further proces^g based on the 
final surfiice profile model and a desired surface profile, to determine process vahies of 
one or more input variables governing a plasma process sequence appropriate for creating 
the desired profile on a process substrata or, by inserting process values of the input 
variables into the final mathematical model, to predictively calculate a process surfice 
profile to be created on a process substrate by a plasma process sequence defined by the 
process values. 

In particular, the hardware system depicted in FIG. 3 may be used to implement 
the calibration procedure illustrated by FIG. 1 as follows. The input variable test values 
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3 selected in step 200, the test value of any fixed input parameter, the rough preliminary 
values of the unknown plasma model coefficients provided in step 215, the rough 
prelinunary values of the unknown profile simulator coeflBcients provided in step 225, and 
the test surface profile created in step 230 and, as needed, the desired surfiace profile 
and/or process values of interest are provided to the database 396 so that they are 

10 available to the analyas module 392, Altemativdy, the module 392 may retrieve any of 
the test values, rough preliminary values and test surface profile data firom the mass 
storage device 357 or user interface 394 in response to a user command. Or, the rough 
preliminary vahies may be detemuned by the module 392, based on the input variable test 
values, according to a predetermined algorithm. 

15 By executing the plasma modding and profile simulation of stq)s 2 1 0 and 220, 

respectively, the module 392 establishes the initial mathematical surfece model pre<Ucting 
the profile created by the test process. In step 240 the module 392 accesses the test 
surface profile and compares it with the initial mathematical sur&ce profile model and 
evaluates the re^dual according to some predetwmined criterioa If the reddual is not 

20 sufficiently small, the analysis module 392 uses the results of the comparison to adjust the 
vahicsofthe unknown plasma modd and profile simulator coeffideiitsm step 250. The 
new values are retained in the database 396 for anothw iteration of the 
modeling/simulation and comparison steps. When the test surfece profile and 
approximate prediction are suffidcntly aniilar, the coeffidwits used in that final iteration 

25 are stored in the database 396 as the optimum values. 

The analysis module uses these optimum values of the input variables for 
computation of process vahies, which can be loaded into a plasma reactor for production 
of a device including the derired profile, or for profile prediction as described above. 

It will therefore be seen that the foregoing represents a highly extensible and 
30 advantageous approach to plasma processing of samconductor devices. The terms and 
expressions employed herein are used as terms of description and not of limitation, and 
there is no intention, in the use of such terms and expressions, of excluding any 
equivalents of the features shown and described or portions thereof; but it is recognized 
that various modifications are possible withm the scope of the invention claimed. For 
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example, the various modules of the invention can be implemented on a general-puq>ose 
computer using appropriate software instruaions, or as hardware circuits, or as mixed 
hardware-software combinations (wherein, for example, plasma modeling and profile 
simulation are performed by dedicated hardware components). 
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CLAIMS 

1. A method of determining a respective process value of at least one input variable 
governing a plasma process sequence for creating a de^ed surface profile on a process 
substrate, the method comprising the steps of: 

a. selecting a respective test value of the at least one input variable; 

b. subjecting a test substrate to a test process defined by the respective test value, 
thereby creating a test sur&ce profile; 

a providing an initial sur&ce profile model in t^ms of the at least one input 
variable and at least one unknown coeffident; 

d. generating an approTdmate profile prediction firom the initial sur&ce profile 
model and the respective test value of the at least one iiqmt variable; 

e. generating an indicator of difference between the test sur&ce profile and the 
jq;)proximate profile prediction; 

£ generating a respective optimum vahie of the at least one unknown coefiBd^ 
that miniimzes the indicator of difference; 

g. modifying the initial sur&ce profile model to indude the at least one optimum 
value, thereby providing a final modd in terms of the at least one input variable; and 

h. generating the respective process value of the at least one input variable firom 
the final modd and the desired sur&ce profile. 

2. The method of claim 1 y^erdn the at least one unknown coeffid^ comprises a 
plurality of unknown coefiBdents. 

3 . The method of daim 1 wherdn the at least one input variable comprises a plurality of 
input variables, the approximate profile prediction being generated fi'om the initial sur&ce 
profile model and the respective test values of each of the plurality of input variables. 

4. The method of claim 1 wherem the test surface profile comprises a plurality of 
snapshots, the approximate profile prediction indudii^ a &Bme corresponding to each 
snapshot, the step of generating an indicator of difference between the test sur&ce profile 
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and the approximate profile prediction including pmrwise comparison of each snapshot 
witKthe respective corresponding frame. 

5. The method of claim 1 wherdn generating an approximate profile prediction includes 
using a respective rough preliminaiy value of the at least one unknown coefficient. 

6. The method of daim 35 wherein generating a respective optimum value of the at least 
one unknown coefficient includes changing at least one of said at least one respective 
rough preliminary value of the at least one unknown coefficient and comparing the test 
surface profile and the approximate profile prediction incorporating the at least one 
changed vahie. 

7. The method of claim 1 wherdn the indicator of diflferaMe is generated by a 
multidimensional nonlinear least-squares technique. 

8. A method of proces^g a process substrate to create a desired surface profile thereon, 
the process compriang the step of applying a plasma process to the process substrate, the 
plasma process being defined by the respective process vahie of the at least one input 
variable determined by the method of daim 1 . 

9. A device made by processing a process substrate as defined in claim 4. 

10. A method of predictivdy calculating a process surface profile to be created on a 
process substrate by a plasma process sequence defined by a respective process value of 
at least one input variable, the method comprising the steps of: 

a. sdecting a respective test value of the at least one input variable, at least one of 

5 said at least one respective test value being unequal to at least one of said at least one 

6 respective process value; 

7 b. subjecting a test substrate to a test process defined by the respective test vahie, 

8 thereby creating a test surface profile; 
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c. providing an initial surface profile model in terms of the at least one input 
variable and at least one unknown coefiBcient; 

d. generating an approximate profile prediction firom the initial surface profile 
model and the respective test value of the at least one input variable; 

c. generating an indicator of diflFwence between the test surface profile and the 
approximate profile prediction; 

£ generating a respective optimum value of the at least one unknown coeflBcient 
that minimizes the indicator of difference; 

g. modii^Jig the initial surface profile model to inchide the at least one optimum 
value, thereby providing a final modd in terms of the at least one input variable; and 

h. introducing the respective process value of the at least one input variable into 
the final model, thereby forming a description of tiie process surfece profile. 

1 1. The method of dann 10 wherein the at least one unknown coeffident comprises a 
plurality of unknown coeffidents. 

12. The method of claim 10 wherein the at least one input variable comprises a plurality 
of input variables, the appro»mate profile prediction being generated irom the initial 
surface profile modd and the respective test values of each of the plurality of input 

4 variables. 

1 13, The method ofclaim 10 wherdntiie test surfece profile conq)Tises a phirality of 

2 snq>shots, the approximate profile description including a fiame corresponding to each 

3 snapshot, the step of generating an indicator of difference between the test surface profile 

4 ai^ tiie approximate profile description including pairwise comparison of each snapshot 

5 with the respective corresponding firame. 

1 14. The method of daim 10 wherein generating an approximate profile description 

2 indudes using a respective rough prdimmaiy value of the at least one unknown 

3 coefiBdent, 
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1 15. The method of claim 14 wherein generating a respective optimum value of the at least 

2 one unknown coeffident includes changing at least one of said at least one respective 

3 rough preliminary value of the at least one unknown coe£Scient and comparing the test 

4 sur&ce profile and the approximate profile prediction incorporating the at least one 
s changed value. 

1 16. The method of claim 10 wherein the indicator of difference is generated by a 

2 multidimensional nonlinear least-squares techiuque. 

1 17. The method of daim 1 0 wherdn the respective process value of the at least one input 

2 variable varies with time. 



1 1 8. A method of configuring an apparatus for processing a process substrate according to 

2 a plasma process sequence defined by a respective process value of at least one input 

3 variable, the apparatus including a plasma reactor, the at least one input variable indudtng 

4 at least one reaction variable, the method conq>rising the steps o£ 

5 a. sdecting a respective test value of the at least one input variable; 

6 b. subjecting a test substrate to a test process defined by the respective test value, 

7 thereby creating a test sur&ce profile; 

s c. providing an initial surface profile model in terms of the at least one input 

9 variable and at least one unknown coefiBdent; 

10 d. generating an {q)proximate profile predictibn firom the initial sur&ce profile 

11 modd and the respective test value of the at least one input variable; 

12 e. generating an indicator of difference between the test surface profile and the 

13 approximate profile pre(tiction; 

14 f generating a respective optimum value of the at least one unknown coeffident 

15 that mimmizes the indicator of difference; 

16 g. modifying the initial sur&ce profile modd to indude the at least one optimum 

17 valuer thereby providing a final modd in terms of the at least one input variable; 

18 h. generating the respective process value of the at least one input variable from 

19 the final model and the desired surface profile; and 
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20 i. configuring the reactor to process the process substrate according to the derived 

31 respective process value of the at least one reaction variable. 

1 19. An apparatus for determining a respective process value of at least one input variable 

2 governing a plasma process sequence for creatmg a desired surface profile on a process 

3 substrate, the apparatus comprising: 

4 a. a compxiter memory for storing the desired surfiice profile; 

5 b. a computer memory for storing a test surfiice profile, created by subjecting a 

$ test substrate to a test process defined by a respective test vahie of the at least one input 

7 variable; 

8 c. means for generating an initial surface profile model in terms of the at least one 

9 input variable and at least one unknown coefficient; 

10 d. means for generating an ^proTomate profile description firom the initial sur&ce 
a profile model and the re^ective test value of the at least one input variable; 

12 e. means for generating an indicator of difference between the test surface profile 

13 and the approximate profile prediction; 

14 £ means for gwierating a respective optimum value of the at least one unknown 

15 coefficient that nmimizcs the indicator of Terence; 

16 g. means for modi^g the initial surfece profile model to include the at least one 

17 optimum value, thereby providing a final model in terms of the at least one input variable; 

18 and 

19 h. means for generating the respective process vahie of the at least one input 

20 variable fi-om the final model and the desired surfiice profile. 

1 20. The apparatus of claim 19 wherm the test surface profile comprises a plurality of 

2 snapshots, the approximate profile prediction indudfaig a fi^ie corresponding to each 

3 snapshot, the means for generating an indicator of difference between the test surface 

4 profile and the approximate profile prediction being configured so as to compare eadi 

5 snapshot with the respective corresponding fiame. 
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1 21 . The apparatus of claim 19 &rther comprising a computer memoiy for storing a 

2 respective rough preliminary value of the at least one unknown coeffident, the means for 

3 generating an approximate profile description from the initial surface profile model and 

4 the respective test vahie of the at least one mput variable employing the respective rough 

5 preliminary value. 

1 22. The method of claim 2 1 wherein the means for generatmg a respective optimum 

2 value of the at least one unknown coefficient is configured to change at least one of said 

3 at least one respective rough prdimmary value of the at least one unknown coeffident and 

4 to compare the test sur&ce profile and the approximate profile prediction incorporating 
s the at least one changed value. 



1 23. The method of daim 19 wherein the means for generating an mdicator of difference 

2 between the test surface profile and the approTdmate profile prediction employs a 

3 multidimensional nonlinear least-squares technique. 

1 24. An apparatus for predictively calculating a process sur&ce profile to be created on a 

2 process substrate by a plasma process sequence defined by a respective process value of 

3 at least one input variable, the method compri^ng the steps of. 

4 a. a computer memory for storing the req[>ective process vahie; 

5 b. a computer memoty for storing a test sui&ce profile, created by subjecting a 

6 test substrate to a test process defined by a respective test value ofthe at least one iiq>ut 

7 variable; 

8 c. means for generating an initial sur&ce profile model in terms of the at least one 

9 input variable and at least one unknown coefficient; 

10 d. means for generating an approTdmate profile prediction from the initial model 

11 and the respective test value of the at least one input variable; 

12 e. means for generating an indicator of ctifference between the test sur&ce profile 

13 and the appro^dmate profile prediction; 

14 f. means for generating a respective optimum value of the at least one unknown 

15 coeffident that minimizes the indicator of difference; 
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g. means for modi^ng the initial surface profUe model to include the at least one 
optimum value, thereby providmg a final model in terms of the at least one input variable; 
and 

h. means for introducing the respective process value of the at least one input 
variable mto the final model, thereby forming a description of the process surface profile. 

25, The apparatus of claim 24 wherein the test surface pirofile comprises a ptarality of 
snapshots, the approximate profile prediction inchiding a prediction corresponding to each 
snapshot, the means for generating an indicator of diflference between the test sur&ce 
profile and the approximate profile prediction being configured so as to compare each 
snapshot with the respective corresponding prediction. 

26. The apparatus of claim 24 fiirttiercompriang a coniputcr memory fi^rstoi^ 

E respective rough preliminary value of the at least one unknown coeflScient, tiie means for 
J genOTting an approximate profile prediction emplojdng the respective rou^ prelinunary 
4 value. 

1 27. The metiiodofclaim 26 wherein Uie means for generating a respective optimum 

2 value of tiie at least one unknown coeflBdent is configured to change at least one of said 

3 at least one respective rough preliminary value of the at least one unknown coefiBd«it and 

4 to compare the test surface profile and the approximate profile prediction incorpora^ 

i the at least one changed value. 

1 28. The mctiiod of daim 24 wherdn the means for generating an mdicator of diflference 

2 between the test surfiice profile and the approximate profile prediction employs a 

3 multidimensional nonlinear least-squares tedmique. 

1 29. The method of daim 1 fiirther comprising the step of applyir« a plasma process to 

2 the process substrate, tiie plasma process being defined by tiie respective process vahie of 

3 the at least one input variable. 
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